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We study constraints on left-right symmetric models from searches of semileptonic decays of D, 

Ds, and B mesons, mediated by heavy neutrinos N with masses niN ~ GeV that go on their 
mass shell leading to a resonant enhancement of the rates. Using these processes we examine, as 
a function of mv and Mw^, the physics reach of the recently proposed high-intensity beam dump 
experiment SHiP, which is expected to produce a large sample of Ds mesons. We compare these 
results with the corresponding reach of neutrinoless double beta decay experiments, as well as like- 
sign dilepton searches with displaced vertices at the LHC. We conclude that the SHiP experiment 
has clear advantages in probing the left-right symmetric models for heavy neutrinos in the GeV 
mass range. 
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I. INTRODUCTION 

One of the current avenues to search for new physics 
has been opened with the discovery of neutrino oscilla¬ 
tions.^ Oscillations indicate that neutrinos necessarily 
must be massive, while the standard model (SM) of elec- 
troweak interactions considers neutrinos as massless par¬ 
ticles. Even if one tries to include the neutrino masses 
within the SM by means of Yukawa interactions —the 
mechanism that gives mass to all other fermions— the 
scenario does not come out natural, because one needs 
to introduce right-handed neutrino fields, which in turn 
lead to the possible inclusion of Majorana mass terms 
for these extra neutrino components, thus breaking the 
lepton number, and so on [2-4]. In summary, the sce¬ 
nario opens into a wide range of possibilities, all of which 
indicate the existence of extra particles, mainly extra 
neutrinos, in a very broad range of energy scales [5].^ 
While each extension of the SM proposes extra neutrinos 
in a specific mass range, in general, the possibilities are 
from a few eV all the way to grand unification scales of 
the order of 10^® GeV. On the other hand, the experi¬ 
ments can put bounds on specific combinations of neu¬ 
trino masses and mixings, each one covering a different 
range of masses [14]. 

In a typical seesaw model of neutrinos based on the 
SM gauge group, the heavy neutrinos are coupled to the 
standard sector through a small mixing with the stan¬ 
dard leptons in the electroweak currents. In contrast, 
in a left-right symmetric gauge theory, the heavy neutri¬ 
nos connect to the standard sector primarily through the 
right-handed currents. In the first case the couplings are 


^ For the current status of oscillation data, see, for example. 
Ref. [1]. 

^ The role of heavy neutrinos in astrophysics and cosmology has 
been studied in Refs. [6—13]. 


suppressed by the mixing, while in the second case they 
are suppressed by the large mass of the Wr bosons. In 
the present paper, we study the latter scenario, applied to 
rare decays of heavy mesons that are induced by massive 
neutrinos that go on mass shell in the intermediate state. 
Depending on the mass of the heavy neutrinos and the 
Wr bosons, experiments could either discover the effect 
or establish new bounds, restricted to neutrino masses 
below the mass of a decaying particle, i.e., below 5 GeV 
for B mesons or below 1.8 GeV for D and Dg mesons, so 
that the resonant effect can occur. 

A particularly attractive feature of left-right symmet¬ 
ric (LR) models of electroweak interactions is that the 
appearance of right-handed neutrino components is not 
accidental but required in order to complete the right- 
handed lepton doublets [15-17]. In these models, the 
heavy neutrinos are mainly the right-handed fields with 
small admixtures of the standard neutrinos. Gonse- 
quently, they couple sizably to the SM sector through 
the gauge fields of the right-handed sector, Wr, but the 
latter induce very weak interactions on the SM particles 
due to the large scale of Mw^ ■ The neutrino mass range 
of interest here is appropriate for searches in B factories 
such as BABAR [18] and BELLE [19] and future high- 
intensity beam experiments such as SHiP [20]. 

This article is organized as follows. In Sec. II we briefly 
summarize the main formulas for production and decay 
of heavy neutrinos in LR models. In Sec. Ill we dis¬ 
cuss the sensitivity of heavy neutrino searches driven by 
meson decays D, Dg, and B in LR models as could be 
observed in the SHiP experiment, comparing them with 
current and future limits coming from neutrinoless dou¬ 
ble beta decay (OuP/S) experiments and equal-sign dilep¬ 
ton searches with displaced vertices at the LHG. We leave 
the last section for a short discussion and summary. 
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II. THEORETICAL SETUP 


Left-right symmetric models are well-motivated and 
popular extensions of the standard model based on the 
gauge group SU (2 )l x SU{2)fixU (1)b-l with gauge cou¬ 
plings gL, Qr, and gi, respectively, and a LR symmetric 
assignment of quarks and leptons: Ql,r = iu,d)L,R and 
Ll,r = {v,1)l,r [15-17]. These theories have two partic¬ 
ularly interesting features: First, LR symmetric models 
can be accommodated in broader groups such as Pati- 
Salam or S'O(IO); and second, these models contain three 
right-handed neutrinos, which are required to complete 
the right-handed lepton doublets, making it a natural 
framework to justify the type-I seesaw mechanism. Thus, 
in a three-family scenario, there are three light (or stan¬ 
dard) and three heavy neutrino mass eigenstates, which 
here we call Ni. 

At some mass scale larger than the scale of 

the electroweak symmetry breaking, the gauge group 
SU{2)l X SU{2)^ X 17(1)b-l is broken down to the SM 
group SU{2)l X C/(1)y. The charged current interactions 
relevant for the present analysis are 

^ = ^ * drpRU+ • irPRN ) 

^ _ _ ( 1 ) 
+ ^ (Kd • Ii^PlN ) , 

where and are the quark mixing matrices, and 
and are the lepton mixing matrices, of the left- 
and right-handed sectors, respectively. In the rest of this 
work, we will assume = 1 for simplicity, and thus we 
will simply denote = U. The charged boson states 
and are linear combinations of mass eigenstates, 
denoted as and W^'- 

Wl = cos C - sin C, 

IF± = IF± sin C -f W:^ cos C, 

where the Wr-Wr mixing angle C, to a good approxima¬ 
tion [21], is given by 


tan 2^ 


‘2gLgRM^^ sin 2(3 
9r^Wl + ~ 


^9r ^Wl 
gL Mwr 


sin 2/3 


(3) 


with tan (3 = k'/k being the ratio of the two vacuum ex¬ 
pectation values of the SU{2) bidoublet Higgs $. Here 
the approximate expression refers to Mwj, ^ ■ The 

masses of Wf 2 are denoted as M\y^ and Mw^, respec¬ 
tively. The maximal value of the Wr-Wr mixing cor¬ 
responds to K = k' when sin 2,5 = 1. The perturbativ- 
ity condition gji/4:Tr < 1 together with gR = g{Mz) ~ 
0.64935 [22] sets the upper bound 


tan2C < 


g{Mz) ■ 


(4) 


The total decay width of each of the three heavy neu¬ 
trinos of the LR models receives a negligible contribution 
from the neutral current interactions and can be written 
in the mass range of uin ^ 1-80 GeV, to a good approx¬ 
imation, as [21, 23, 24] 



Here the masses of all final state particles are neglected. 
For completeness we give the expression for the half-life 
Ti /2 of a neutrinoless double beta {0i'/3l3) decay mediated 
by Wr and heavy N exchange [21, 23, 25]: 


r-/i =Goi \Mn? 


mp 

gl 


( 6 ) 


where Gqi is the phase space factor and AIat is the nu¬ 
clear matrix element. 


III. PHENOMENOLOGY OF LR SYMMETRIC 
MODELS 


For simplicity, we will first restrict our discussion to 
heavy neutrino mixing with the electron sector only. We 
assume that the production of a heavy Majorana neu¬ 
trino in a meson decay, M —)■ eN (with M either D, Dg, 
or B), and its subsequent decay as N ^ eir, is dominated 
by the Wr boson exchange, as is shown in Fig. 1. 



FIG. 1. Heavy neutrino production (a) and decay (b) in the 
LR symmetric model. 


Consider a beam of heavy neutrinos N with average 
energy G, produced in a fixed target experiment. In the 
manifestly LR symmetric scenario, i.e., gR = gR, the de¬ 
cay length of such neutrinos can be written according to 
Eq. (5) as a function of the masses ttin and mwR [25] in 
the form 


L = c-/T^ 


/IGeV 
V rn-N 


/ mwg Y 
\neYj 


(7) 


with 7 = E/mN- Therefore, for mwn ~ TeV and itin ~ 
GeV this decay length is rather large, and one can expect 
that only a fraction of the heavy neutrinos will decay 
inside a given detector. 

In a LR symmetric model, the number of signal events 
in the form of meson decays M —>■ eN followed by the 
heavy neutrino decay V —>■ ctt are suppressed by a factor 
of {mwL■ This suppression comes from the pro¬ 
duction of heavy neutrinos M —> eN, which is suppressed 
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by (rawL !'ynw 1 and from the fraction of heavy neutri¬ 
nos decaying inside the detector, which is proportional 
to also suppressed by a factor of (tovVi,/”^Wh)^- 
From the experimental point of view, this signal cannot 
distinguish a LR symmetric model from a seesaw model 
based on the SM gauge group, where the number of signal 
events is suppressed by the small heavy-to-light neutrino 
mixing |C/eArp x |C/eAfp [20, 24, 26-28] (see also [29-38]). 
However, for the same reason, if limits on |17ew|^ are 
known experimentally, we can extract limits on mwa , by 
doing the conversion: 







( 8 ) 


in the case of the decay Dg —t eN —)■ ecTr. For the 
B —> eN —>• ecTT one should replace Vcs by Vub in the ex¬ 
pression above. It is expected that the quark mixing ma¬ 
trices are similar ~ [39, 40]. Furthermore, in 

the manifestly LR scenario, i.e., gn = g^, the conversion 
reduces to l?7eAr|^ —t {’mw^/innwR)'^■ 

The BELLE experiment, using B meson decays B 
XIN followed by t CTT, has set limits on \UeN\’^ 
10“^ in the heavy neutrino mass range 0.5 —5 GeV [19].^ 
In addition, the proposed high-intensity beam dump ex¬ 
periment SHiP [20], which should produce a very large 
number of Dg mesons, will be sensitive to \UeN\’^ 
10“® — 10“^° for heavy neutrino masses of the order 
of TOjv I GeV [14, 26, 27]. In the same way, these 
BELLE searches on heavy neutrinos correspond to lower 
limits on itiwr near 800 GeV, which are not better than 
the current limit mwR > 2.5 TeV coming from Kq — Kq 
mixing [23, 39, 40, 42-44]. On the other hand, the fu¬ 
ture SHiP experiments will be sensitive to larger 
masses, up to mwR 8—18TeV, for heavy neutrino 
masses m^r ^ 1 —2 GeV, as shown in Fig. 2. 

The LHG can also constrain LR symmetric models 
by searching for like-sign leptons plus two jets [45-48].^ 
Gurrently, the LHG has imposed a lower limit for mwR 
at 3.0 TeV [48] in a neutrino mass range uin ~0.2- 
2.0 TeV. These limits are expected to be extended up 
to mwR ^ 6 TeV [56] in future searches. Although 
the current searches at the LHG are not sensitive to 
heavy neutrinos with masses mjv as low as a few GeV, 
it is still possible to search for heavy neutrinos with 
masses ttin < 80 GeV using displaced vertices [25, 57- 
60].® These heavy neutrinos could be produced via a 
Wn boson as is shown in Fig. 3. The advantage of a 


® The LHCb experiment has searched for heavy neutrinos by using 
the B meson decay mode B —>■ fiN followed hy N fin in a mass 
range mjy ^ 0.5-5GeV, setting limits on |r/^ivP ~ 10“® [41], 

^ Generic processes of this kind have been discussed in Refs. [49— 
52] (see also [53—55]) as tree-level high-energy completions of lep¬ 
ton number violation operators that generate neutrinoless double 
beta decay. 

® Recently, the authors of Refs. [61, 62] have proposed to search 
for heavy neutrinos induced by Higgs decays at the LHC. 




FIG. 2. Regions in parameter space, which can be probed 
by a displaced vertex search at the LHC, QvPP, and SHiP 
experiment, (a) Top: mwR vs niN for fixed qr = Ql- (b) 
Bottom: QrIql vs mwR for fixed miy = 1 GeV. For details, 
see the text. 


displaced vertex search is that, for a vertex separation 
roughly between 10“® and 10° m, there is little or no 
background [25]. A signal based on this displaced vertex 
in principle could be affected by backgrounds from rare 
displaced hadron decays or from pileup effects [58, 63]. 
These backgrounds are found to be negligible for high 
track multiplicity and high enough px cuts, which in our 
case do not eliminate the signal due to the large mass of 
the Wr. 

Using the narrow width approximation for the inter¬ 
mediate N propagator, we can express the number of 
events of this type at the LHC decaying in the range 
from di = 1 mm to d 2 = 1 m as 

X a{pp eN) x Br{N — >■ ejj) x Pjy. (9) 

Here Pjy = [exp(—di/L) — exp(—(i 2 /L)] is the probability 
for a heavy neutrino to decay within a distance between 
di = 1 mm and d 2 = 1 Ta from its production point. 
Following the analysis of Ref. [25], we have used Eq. (9) 
to estimate the parametrical region that would produce 
at least five events with vertex separation between 1 mm 
and 1 m, for an integrated luminosity of £ = 300 fb“^ and 
collision energy y/s= 13 TeV, and using the kinematical 
cuts in transverse momentum and rapidity pt> 40 GeV 
and g < 2.5, respectively, for both leptons and jets. 

In Fig. 2(a), we show two limits from the nonobser¬ 
vation of Ovjjjd, assuming that the contribution via the 
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FIG. 3. Top: Quark-level Feynman diagrams for left-right 
realization of neutrinoless double beta decay with heavy neu¬ 
trino exchange. Bottom: The corresponding diagram at the 
LHC. 

Wn boson and heavy neutrino exchange shown in Fig. 3 
dominates. The weaker limit (solid line) [23, 64-66] cor¬ 
responds to the current bounds on the half-life of neu¬ 
trinoless double beta decay, Ti/ 2 > 2 x 10^^yr [67, 68], 
while the stronger limit (dashed line) corresponds to an 
expected future sensitivity of ri/ 2 > 10^^yr. The grey 
band bordered by the solid line shows the region in pa¬ 
rameter space where a displaced vertex search at the LHC 
could yield at least five events. The solid line near the 
bottom of Fig. 2(a) shows the region in the parameter 
space which can be probed by the heavy neutrino search 
at SHiP. This is so far an unexplored domain of the 
niN — mwn parameter space below toat ~ 5 GeV where, 
as seen from Fig. 2(a), the SHiP searches are much more 
sensitive than the LHC and even more sensitive than op¬ 
timistic future 0iz/?/3 results for heavy neutrino masses 
rriN ~ 1 —2 GeV. 

Up to here, we have made the assumption that ga = 
gL, which corresponds to the manifestly left-right sym¬ 
metric scenario, a case that may not necessarily be so in 
general [69]. In Fig. 2(b), we show, just as in Fig. 2(a), a 
comparison between the Ojz/ 3/3, LHC, and SHiP sensitiv¬ 
ities, but for a fixed heavy neutrino mass rriN = I GeV 
and in the plane gn/gh — ttiwr- As shown, and in agree¬ 
ment with the previous analysis, the LHC is not compet¬ 
itive with Oz//3/3 or SHiP for itin = 1 GeV. It is also clear 
that the SHiP experiment is capable to probe a much 
larger area of this parameter space than even future neu¬ 
trinoless double beta decay experiments. 

The lepton flavor violating muon decays fi —>■ ey, 
/i —> eee, and /r —>■ e conversion in nuclei® also im¬ 


® The authors of Ref. [70] have studied the contribution of Majo- 


pose constraints on LR symmetric models. However, for 
masses itin ^ GeV these bounds are very weak and not 
comparable to those coming from Ojz/3/3 [75, 76]. 

Finally, we close by mentioning that so far we have con¬ 
sidered heavy neutrino mixing only in the electron sector. 
For heavy neutrino mixing in the /t sector, the LHC and 
SHiP limits shown in Fig. 2 will remain the same, but 
the limits coming from Qvj3j5 are not applicable. This is 
so because the OvjSP amplitude is sensitive only to the 
UeN mixing matrix elements. For heavy neutrino mixing 
predominantly in the r-lepton sector, the situation is also 
different. On the one hand, the limits from the LHC will 
be much worse than the limits shown in Fig. 2, due to the 
relatively poorer r reconstruction efficiencies at the LHC. 
On the other hand, the SHiP experiment could be sen¬ 
sitive to I Ur AT I ^ ~ for heavy neutrino masses 

of the order of to at 1 GeV [14], which correspond to 
strong limits on the Wr boson mass up to mwR ~ 9 TeV 
for heavy neutrino masses to at 1 GeV and gn = gL- 

IV. SUMMARY AND CONCLUSIONS 

We have analyzed the scenario of heavy neutrinos 
within left-right symmetric models of electroweak inter¬ 
actions. We considered both a manifestly left-right sym¬ 
metric case, as well as cases where the right and left 
gauge couplings are different. Concerning the heavy neu¬ 
trino sector, we studied the cases where one neutrino is 
in the mass range toat from 1 to 80 GeV, and the rest 
are assumed to be heavier or not to interfere with the 
processes we consider. In LR symmetric models of neu¬ 
trinos, the heavy neutrinos couple to the standard sec¬ 
tor mainly through the right-handed currents, so their 
suppression is due to the large mass of the Wr bosons. 
We studied the sensitivity of Ojz/ 3/3, LHC, and SHiP ex¬ 
periments to the mass of the heavy neutrinos and Wr 
bosons. In the case of LHC, we considered signals of 
equal-sign dileptons with displaced vertices, which are 
free of backgrounds, and in the case of SHiP, we consid¬ 
ered the production of heavy neutrinos through charmed 
meson decays. We find that in the range of masses where 
SHiP is relevant, its sensitivity could be much stronger 
than the other two kinds of experiments, with the poten¬ 
tial of discovery or considerable improvements on current 
niwu bounds. As a final comment, we should 

add that it could be possible to distinguish at SHIP a 
type-I seesaw scenario from a LR symmetric scenario by 
comparing the rates of semileptonic (i.e., N —)■ ctt) with 
purely leptonic (i.e. N —)• eev) modes. Because of the 
strong suppression of v production through the right- 
handed current, in the LR symmetric scenario the purely 
leptonic mode is highly suppressed, while in the type-I 
case the two modes are comparable. 


rana neutrinos in /r —>■ e conversion in nuclei. See also Refs. [71— 
74]. 
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